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ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ 


УДК 539.3 


Определение изменения формы поверхности непрерывно-неоднородного 
термоупругого полупространства при локальном нагреве! 


Л. И. Кренёв, С. М. Айзикович, Б. И. Митрин 
(Донской государственный технический университет) 


Рассматривается осесимметричная квазистатическая задача термоупругости для Ффункционально- 
градиентного полупространства, модуль упругости, коэффициент Пуассона, коэффициенты теплопроводно- 
сти и линейного расширения которого непрерывно изменяются в приповерхностном слое. Предполагается, 
что область внутри круга нагревается источником тепла с постоянной во времени температурой. Вне круга 
поверхность идеально теплоизолирована. Для решения задачи используются аналитические методы, в част- 
НОСТИ, аппарат интегральных преобразований Ханкеля. Решение задачи на первом этапе сводится к реше- 
нию краевой двухточечной задачи для системы обыкновенных дифференциальных уравнений с переменны- 
ми коэффициентами шестого порядка. Для организации устойчивого процесса численного построения реше- 
ния системы обыкновенных дифференциальных уравнений используется метод модулирующих функций. 
В результате решение смешанной граничной задачи сводится к построению решения парного интегрального 
уравнения. Свойства трансформанты ядра интегрального уравнения задачи позволяют применить хорошо 
обоснованный и развиваемый в настоящее время двусторонний асимптотический метод. С помощью данного 
метода найдены в аналитическом виде приближённые выражения величины теплового потока и смещения 
поверхности полупространства. Приведены численные результаты, отражающие искривление поверхности 
неоднородного полупространства для различных случаев изменения механических и температурных свойств 
в приповерхностном слое под действием равномерной температуры в пределах единичного круга. Рассмат- 
риваются случаи, когда значения термоупругих свойств покрытия совпадают со значениями термоупругих 
свойств подложки, либо когда значение характеристики отличается в 2 раза (в большую или в меньшую сто- 
рону) на поверхности и линейно убывает (или растёт) по глубине до значения характеристики в подложке. 
Показано, что максимальное влияние на величину максимального выпора поверхности оказывает разнона- 
правленное изменение коэффициентов теплопроводности и линейного расширения в покрытии. 

Ключевые слова: смешанные задачи, неоднородные материалы, термоупругость, функционально- 
градиентные материалы, аналитические методы 


Введение. Учёт неоднородности свойств материала при моделировании тепломеханического 
воздействия является актуальной задачей термоупругости. Такая задача возникает, например, 
при исследовании свойств защитных покрытий из функционально-градиентных материалов для 
частей машин и механизмов, подвергающихся интенсивному тепломеханическому воздействию. 
Подробное изложение основ термоупругости представлено в монографиях В. Новацкого [1] и 
А. Д. Коваленко [2]. Среди тех, кто активно занимается решением контактных термоупругих задач 
для однородных тел, стоит отметить Дж. Р. Барбера. В его работе [3] приводятся аналитические и 
численные методы решения задачи термоупругой устойчивости, а также приводятся результаты 
по изучению термоупругих деформаций. 

В последнее время в отечественной и зарубежной научной литературе наблюдается рас- 
тущий интерес к решению задач термоупругости для неоднородных материалов. Для задачи теп- 


' Работа выполнена при поддержке Министерства образования и науки РФ (ГК № 11.519.11.3028, соглашение 


№ 14.В37.21.1632), а также грантов РФФИ (12-07-00639-а, 13-08-01435-а). 
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лопроводности решения для отдельных видов законов неоднородности (экспоненциальный, ли- 
нейный законы) приведены в работе [4]. Задача для слоистого покрытия решается в работе 
И. Г. Величко, И. Г. Ткаченко [5]. Собственный метод решения смешанных задач для экспоненци- 
ального закона неоднородности покрытия развивает М. Мода [6—7]. Решение плоской задачи для 
термоупругого взаимодействия в присутствии трения рассматривали также Ци, Ке, Мапа [8]. 

Аналитические методы для решения термоупругих задач применяет Краснюк П. П. Так, в 
работе [9] рассмотрено решение плоской задачи контактного взаимодействия жёсткого теплопро- 
водного штампа кругового поперечного сечения с упругим слоем. 

Решения, используемые в приведённых выше работах, не могут быть применены для 
сложных случаев изменения термоупругих и механических свойств по глубине слоя. 

В настоящей работе рассматривается осесимметричная квазистатическая задача термо- 
упругости для функционально-градиентного полупространства, модуль упругости и коэффициен- 
ты Пуассона, теплопроводности и линейного расширения которого непрерывно изменяются в 
приповерхностном слое. Для её решения используется двусторонне-асимптотический метод, из- 
ложенный в статье [10]. 

Постановка задачи о воздействии локального источника с постоянной температурой 
на непрерывно-неоднородное полупространство. Рассмотрим полупространство @, термо- 
механические характеристики которого непрерывно меняются с глубиной в пределах прилегаю- 
щего к поверхности слоя толщины Я, а затем стабилизируются и остаются постоянными. С полу- 
пространством свяжем цилиндрическую систему координат (г,ф,2). При этом полагаем, что 


О<г<о, а 2<0. Обозначим через и,и,и смещения вдоль осей к ФФ, д а через 


о, 0, О, То’ Г — радиальное, угловое, нормальное и тангенциальные напряжения соответ- 


гф! 7! Тя 


ственно. Кроме того, введём обозначения: для температуры — 7’, коэффициента теплопроводно- 
сти — Л, (2), коэффициента линейного расширения — а, (2), коэффициента теплоёмкости — 
ст(2). В ненапряжённом состоянии температуру полупространства примем равной 7, и будем 
анализировать разность температур 9 =Т -Т.. 


На площадке 0 <г<а поддерживается постоянная температура и образуется тепловой 
поток, направленный внутрь неоднородного полупространства. 
Далее будем полагать, что коэффициенты Ламе М(2) и ^(2), коэффициенты теплопро- 


водности ^, (2) и линейного расширения а; (2) являются непрерывными функциями координаты 


7 ‚ такими, что 


(1) 


2е(-©,0) 
( 


тм М 
2«(-о,0) 


шлЛ 
И 


(2)> Ра 


тт @: (2) >а,. >0, тах а; 


2(-°,0) 


) 
ши Л(2)>Л. >0, тах Л(2 
( 


2)>М. >0, тах М(7 


2е(-ю,0) 
2е(-ю,0) 
> 0, пах: ( 
( 


2Е(-—®,0) 
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где Н — толщина неоднородного слоя, сцепленного с подстилающим полупространством, то есть 
глубина, с которой мы полагаем термомеханические характеристики полупространства посто- 
янными; индекс 5 соответствует подстилающему однородному полупространству; С — не- 


однородному слою; Л,, Л", М,, М, Л., Л;, а., а. — произвольные константы. 
Т Т 


Наряду с парой коэффициентов Ламе для описания упругого поведения твёрдого изотроп- 
ного тела используются: модуль сдвига С и коэффициент Пуассона у или модуль Юнга Е и ко- 
эффициент Пуассона у. Коэффициенты Ламе Л и М (иногда обозначаемый Си называемый мо- 
дулем сдвига) связаны с модулем Юнга Е, коэффициентом Пуассона у и изотермическим моду- 
лем объёмного расширения К’ соотношениями: 

вы Е _ Еу _ М(2мМ+ Л) ты Л 
2(1+у)° (1+%)(1-2)'’ Мл о 2(м+№' (2) 
ЗК (2) =3^(2)+2М(2). 
Граничные условия, при сделанных предположениях, имеют вид: 
9=60,, г<а, 
о о. ) 
07 
На границе сцепления неоднородного слоя с однородным полупространством, при 2 =-Н , в силу 
непрерывности, должны выполняться условия сопряжения по смещениям, их производным, тем- 
пературе и тепловому потоку: 
0% (г,-Н)= о (г,-Н), те (г,-Н) =те (г,-Н), 
иб (",-Н) = и? (г,-Н) кие (г,-Н) =и? (г,-Н), (4) 


65 (,-Н) =6° (^,-Н), (8) (^,-Н) = (65) (=,-Н). 
На бесконечности при (г,-Н) — о смещения, деформации и напряжения исчезают. Значения 


разности температур и теплового потока при этом также стремятся к нулю: 


ИВО тг 0:0") = 0, (5) 
ИИ (М, Е, ЕЕ, „› 0,0, 0„,Т,,6,0') =0. 
Требуется определить смещение поверхности покрытия: 
09 
—| =В(г),г<а. (6) 
02 2=0 


А 














Рис. 1. Постановка задачи об изменении формы поверхности непрерывно-неоднородного термоупругого полупространства 
при локальном нагреве 
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Построение решения несмешанной задачи о термическом воздействии на непрерыв- 
но-неоднородное полупространство. При сделанных предположениях соответствующая кра- 
евая задача статической осесимметричной термоупругости в форме Ламе имеет следующий вид 
(без учёта кручения): 


мати (ме ем к), 
м(г) ум + (М (2) + (2) Ем (2) 5+ (= К (26+ (2) ©) 


А; (2)У?8 + Л, (2)0'=0. 


Здесь используются следующие обозначения: 
Аам(2 амл(27 2 
М'(2)= а Ру м А — 
а [вр гог\ дг 07 
ал. (7 ак (2 
О И в ди А [Е )= а. } (=) ( ) 
д г 02’ [@р4 
Будем искать решение для смещений и,и/ и температуры 0 в виде интегралов Ханкеля: 

















ог.) = [И (ь,) (ума (2) = М (2) (меч, (8 
6.2) = [7 (9,2) 2 (ук)чау 


Подставим (8) в систему дифференциальных уравнений в частных производных (7) и, приравняв 
к нулю подынтегральные выражения, получим систему обыкновенных дифференциальных урав- 
нений второго порядка: 
МИ" +у(М+л)И/' - у? (2М ^)и Е МИ’ + УМ'И/ = АУТ, 
(2м+л)И/” - у(М+ Л)’ - уми/ + (2м' + ^)И/'- УЛ’ =КТ +АТ’, (9) 
А (т°) эта] (1) 50, 
где ' указывает на дифференцирование по 2. 
Граничные условия (3) примут следующий вид (для общности мы учитываем также воз- 
можность приложения в пределах круга радиуса а горизонтального воздействия 4 (г) ). 
в рн 
)- [в В(р) 7, (ру)рар. 





(10) 


Используя векторное представление для трансформант: 
Т 
аа (11) 
х, =О, Хх, =0',Х. нии, Х. =И/',Х,; =Т, Хь ЕТ’, 
запишем систему (9) в матричном виде, при этом явно выделим части соответствующие покрытию 
и подложке: 





С 

О = АС, -Н<7<0, (12) 
5 

м = А°х°, -®<2<-Н, (13) 
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Граничные условия при этом имеют вид: 
хв (\,0) =-В(У), 
х. (\,0) = хг (у,-Н), х з (\;0)= хз (у,-Н), 
хё (у,0) =х$ (\,-Н), хё (у,0) = к (14) 
хз (у, 0) = хз (у,-Н), хь (у,0) = хь (у,-Н). 

Необходимо отметить, что последнее условие имеет место в силу непрерывности изменения 
свойств слоя и основания (4). 

Общее решение системы (9) для однородного полупространства Л’=М' =К'=^Л! =0, 


[ет 





М> 0, Л> 0, К >0, Л, >0 имеет вид: 








еее 
хз (у,2) = (а, +(1+у2)9, г к, +к,)2) д, е*, 
У. (У,2) = (4, + (-к, +2) а, + (-к,/у- (к, +к,)2)а а (15) 


хз (у,2) = (4, + (1-к, +у2)а, + (к ‚№ на ее 
х$ (у,2) =а.е*, хё (у,2) = убе*, 
Л+3М К _ К 


Ам" 2 4(2М+Л) 3 4м' 


где 4, (/=1,2,3) — произвольная функция параметра у. 





Решение хб (у,2) системы дифференциальных уравнений (12) строится методом модули- 


рующих функций. Будем искать х° (у,2) в виде: 











Е] 
х° (у,2) = У. а, (\)а, (у,2)е*. (16) 
1=1 
Векторы а, (у,2), (/=1,2,3) определяются из решения следующей задачи Коши: 
а, дса -уа,-Н<2<0,/=1,2,3, (17) 
[64 
при начальных условиях, для 2 =-Н: 
а, (у,2) 2=н = (1, у, 1, у, 0,0) ' 
а, (У,2) 2—н = (У2, У Е у?2, с? К, + уг, У — к, у + у2,0, 0). Г 


а, (у,2 |, м = (-Кь/\ + (к + кз) 2, + (К +к)у2, -Кь/\ + (к + кз) 2, Кь + (к» +к, У?) 


-н" 


Константы а, (у) (7 =1,2,3) определяются из условия (14). 


Таким образом, мы имеем 


2.4, (10 (у) =-В(У), 
м, (у) =м(0о)а, (у,0) +М(0)уа; (у,0), (18) 
м (\) = гоу (+,0) + (^() +2м(@)а' (4,0) & (6)а? (5,0), 

О, (у) =а? (\,0), О, (у) = 0, (у) =0, 


где ах (у,2), (/=1,2,3) обозначает А-ю компоненту вектора а, (\у,2), /=1,2,3, К =1,2,3,4,5,6. 
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о, -® (8, 9 
4; =-КЗВ (У). 
3) _ Аз (У) _1 
ми А,› (У) О (У)' 
зил _ Ав(У) 1 
и А,› (У) О (У) 
& (у) = тт А (У) = М (Ум, (+) - м, (УМ (У) 
О, (У) 
Окончательно получаем следующее выражение для компонент вектора решения х° (у,2), при 
>-Н: 
х, = (у,2) = (В (№2) В(\))е" |, 
х, =" (у,2) = (В (у,2)В (у) е* м, 
х; =И/ (у,2) = (В (у,2)В (у))е" /\ 
х. =И/"(\,2) = ( . (у,2)В (у)}е* М, (20) 
тн.) (в (де 
х; =Т'(у,2)=Е (у,2)В (у)е" /\, 





Введём обозначения: 
Б (",2)= [в (у,2)В (у)е* 1, (уг)ау, (К =1,2 1=1),(К =3,4,5,6 /=0), 
а (21) 
25 (г,2) = [6 (у,2)В (у) е“ 2, (уг)уау, (К =21=0),(К =4 7=1). 
0 


В соответствии с (8), (16) мы можем теперь выписать выражения для смещений, температуры, 
теплового потока и деформаций. 
1 3 
пера ЕЕ (2), 


00(г,2) -18(.2). (22) 





9(г,2)=Ть (г,2), 


Построение замкнутого приближённого численно-аналитического решения смешан- 
ной задачи статической теплопроводности о нагреве неоднородного полупространства 
источником с постоянной температурой в пределах круга радиуса а. Рассмотрим задачу о 
нагревании функционально-градиентного полупространства с поверхности источником с постоян- 
ной температурой в пределах круга радиуса а. Из (21)—(22) следует, что 


09(г,2) = [6 (у,2)В (у) 7 (уг)е*ау. (23) 
На основании (3), (23) получаем, при >. 0: 
[5 (+,0)в(ч)л (уг)е“ау=6,, г<а, (24) 
где 9, — температура т. тепла, а 


В(г) = [в (а)-7. (га)ада, В (у) = [В (р). (Ур) рар. (25) 
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Сделаем замену переменных и обозначим: 


УН =и; Л=Н/а; г’=г/а; р’=р/а; В’ (у)=В (у,0), В'(р’а)=В(р), (26) 
далее штрихи опускаем. Получим, что поставленная задача сводится к решению интегрального 
уравнения вида: 





1 со 
[В(р)рар[ Е (ил) (игл" ) № (ир^” аи =6,, г <1. (27) 
0 0 
Выражение (27) представим в виде парного интегрального уравнения: 
м Л (уг)е" ау =6, г<1, 
(28) 
Г В(у о (ГУ) )уау= 0, г>1. 
Используя операторы 
а‘ ФГ). | : 
= т 
= ь = 3» (г\) = со$Ёу; И!1=1; (29) 
ы а 
= [- Е. 1, (гу) = У" созёу, (30) 
ог - в’ 
преобразуем (28) к следующему виду: 
Тв(у) 18 (Лу) соз (у) Чу =1, Ё<1, 
(31) 


Гв(у) ) сов (у)ау =0, #>1. 


Аналогично работе [11] можно исследовать асимптотические свойства функции /? (и) и показать, 
что её можно аппроксимировать выражением вида [12—13]: 


(и) = В (и)+Е (и). 


Здесь обозначено 








Рассмотрим простейший случай, когда [2 (и) имеет вид: 
2\-2 
У. + А”А- 
Б (№) = В2Л- 2" 
Тогда уравнение (31) для случая, когда область внутри круга некоторого радиуса нагревается 
источником тепла с постоянной по времени температурой, можно записать в виде: 





2-2 


Е у? + А?Л 
В(у)-————— -со$Ёуау =1, Ё<1, 
и о 





- (32) 
[В (у)созёуау =0, #>1. 
Введём обозначение 
= [в (У) созёуау. 
Из уравнения (32) с учётом чётности И части поли 
В(Ё) = СА В (33) 


А2Л-2 ` 
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Используя (33), (32) можно переписать: 





со В? 
[В (у) соз мау = С сп АЛ + Рх, ЕЕ 
: (34) 
[В(У)соз Иду =0, Е>1. 
0 
Обращая преобразование Фурье, имеем: 
2\-2 с; | 1 
ВЕ” ие ОИ ЕОАА + АХ" с0$ и $1 АА” (35) 


п Дл? и п и? + А*Л? 
Постоянную С определим из условия, что В(и) должно удовлетворять исходному уравнению 
(13). Получим: 
ВА" сн АЛ" + АА" $ АА" ВА" 








С В2\-2 — А? ыы А? =0. (36) 
Формулу для теплового потока (25) приведём к более удобному виду, используя равенство Пар- 
севаля: 
111 1 сп АА" ЕН АЛ ЧЕ 
В(г) = сд ме (37) 
(Г) п |2 (0) л-^? М - г? УЕ? г? 


Определение формы изменения поверхности неоднородного полупространства, 
нагретого источником с постоянной температурой в пределах круга радиуса а. Проана- 
лизируем деформацию поверхности неоднородного полупространства под действием локального 
статического нагрева. Согласно (23) при отсутствии поверхностного механического воздействия 
получаем: 


В [В(р)рар [в (\,2).7. (Ур) 7. (уг)е“ ау. (38) 


В случае однородного полупространства формула (38) позволяет записать простые аналитические 
выражения для описания вертикального смещения поверхности полупространства при действии 
постоянного теплового потока и постоянной температуры. Рассматривается установившийся про- 
цесс подвода тепла к полупространству при отсутствии его отвода через поверхность, поэтому 
форма поверхности описывается логарифмической функцией и, для определённости, будем фик- 
сировать некоторое значение г, на оси г, при котором и’, (^,0) =0. При отображении получен- 


ных значений для перемещений поверхности на рисунке 2 примем г, = 3,6. 


1 
ив (г) =и? (г,0)-и? (^,0) т Ее 
т 2 





‚ (Ур) рар Г уг) - 2% (ук) ау = 


[мг -№(и%) )зпуау ыы -1+1- 2 (9%) эпуау = (39) 


2 


ий и, <1, 
мо, ПТ, 
Г 
На рисунке 2 представлено относительное изменение формы поверхности неоднородного полу- 
пространства под действием постоянной температуры, приложенной в пределах круга единичного 
радиуса. Обозначения кривых формируются следующим образом. Последовательно описывается 
характер изменения модуля Юнга, коэффициента Пуассона, коэффициента теплопроводности и 


коэффициента линейного расширения в покрытии: 
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0 — значение характеристики в покрытии не отличаются от соответствующих значений в 
подложке; 

3 — значение характеристики на поверхности в 2 раза превышает значение в подложке и 
линейно убывает по всей толщине покрытия до значения в подложке; 

4 — значение характеристики на поверхности в 2 раза меньше значения в подложке и ли- 
нейно возрастает по всей толщине покрытия до значения в подложке. 

Таким образом, обозначение 0_0_0_0 соответствует однородному термоупругому полупро- 
странству. 


0,030 





0,025 


0,020 





0,015 


и(г)-и(го) 


0,010 


0,005 





0,0 0,5 1,0 1,5 2,0 2,5 3,0 Э,5 4,0 
Г 


Рис. 2. Термоупругое искривление поверхности неоднородного полупространства для различных случаев изменения 
свойств в приповерхностном слое под действием равномерной температуры в пределах единичного круга. Толщина слоя 
равна радиусу круга 


Анализ графиков на рисунке 2 показывает, что наибольшее влияние на величину макси- 
мального выпора поверхности оказывает разнонаправленное изменение коэффициентов тепло- 
проводности и линейного расширения в покрытии. 

Заключение. С помощью двусторонне-асимптотического метода построено решение задачи о 
локальном нагреве поверхности непрерывно неоднородного термоупругого пространство. Опре- 
делены численные значения выпора поверхности при различном характере изменения упругих и 
термомеханических характеристик приповерхностного слоя. Показано, что характер неоднородно- 
сти приповерхностного слоя оказывает существенное влияние на величину выпора поверхности. 
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ОРЕТЕКМТМАТТОМ ОЕ 5УВРАСЕ ОЕРОВМАТТОМ ОЕ СОМТТМООЦ$ЕУ 
ТМНОМОСЕМЕОЦ$ ТНЕКМОЕГА$ТТС НАЕЕ-5РАСЕ ИМОЕВ ОСА НЕАТ!МС' 


Е. Г. Кгепет, $. М. АхНомсн, В. Т. МИгт 
(оп Зае Тесйтса! УпМ№ег$Ку) 


Тре ажа/у зуттейтс диа айс ргоМет оп {Пегтое/азНсйу Юг {те ипсвопа/у дгадеа ва! -зрасе, иЛЁН е/азбсйу 
тоадии$, РоГ5оп гаво, пеаЁ сопдисвоп, апа пеаг ехрапхюоп соетаетё уагутда сопёптиоизу т пе пеаг-зиГасе 
[ауег [5 сопу!аегеа. ГЕ 15 5ирро5еа {таЁ те агеа п те сгсе ийт те Кпоит гад!из 15 Бетд Весеа Бу те ВеаЁ 
5оигсе ий" Ите-сопапЕ {етрегайиге. Ошё5усае те стае, те зипасе 15 регесНу тзиееа. Апау#с тетосб, раг- 
Исш!апу, Напке! тёедга! {гап$Югт, аге изеа 0 зо№е {те ргоШет. АЁ Пг5Е, те 5о/ивоп /5 гедисед Ю о№Мтд а ио- 
ротЁ Боипаагу рго Мет юг те ОБЕб зубет ий" уагутда соейс!епёЕ5 о! те 1 огаЕГ. ТВе тоашайпа Гипсйоп5 
тепод [5 изеа ю ргомае а а Ме питегса! зо/ивоп © {те ог4тагу @егепва! зу$ет. Аз а гези! те оийоп юЮ 
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